On the stability of precessing superfluid neutron stars 
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We discuss a new superfluid instability occuring in the interior of mature neutron stars with 
implications for freely precessing neutron stars. This short-wavelength instabihty is similar to the 
instability which is responsible for the formation of turbulence in superfluid Helium. Its existence 
raises serious questions about our understanding of neutron star precession and complicates attempts 
to constrain neutron star interiors using such observations. 



Introduction. — Neutron stars tend to be extremely sta- 
ble rotators, with stability that sometimes rivals that of 
the best atomic clocks. Yet a growing sample of pulsars 
exhibit spin irregularities, like glitches and timing noise. 
They may also be undergoing free precession. From the 
theory point of view, one might expect precession to be 
generic. Nevertheless, for reasons still to be understood, 
compelling evidence for long-period precession has only 
been found in the timing data of a few pulsars. The best 
candidate is PSR B1828-11 [l] which exhibits a - 500 d 
high-quality periodicity, with an amplitude of a few de- 
grees. The paucity of precessing neutron stars is one of 
the mysteries of pulsar physics. To explain why preces- 
sion is so rare is difficult. After all, a description of even 
a modestly realistic neutron star requires the fusion of 
much of modern theoretical physics. One would need 
to account for strong gravity, supranuclear density mat- 
ter, superfluidity/superconductivity and potentially very 
strong magnetic fields. 

In the standard picture of a mature neutron star the 
bulk of the neutrons are superfluid and rotate by forming 
a dense array of vortices. Meanwhile the outer core pro- 
tons are expected to form a type II superconductor, with 
the magnetic flux carried by fluxtubes. The coupling 
between these two distinct fluid components is usually 
assumed to have the same form as in the case of super- 
fluid Helium, see ^^^l for a recent discussion. However, 
this model is based on the assumption that the neutron 
vortex array is (locally) straight. This may not be the 
case. In a body that undergoes a more complex motion 
one might expect to find that the vortices get tangled 
up to form a state of superfluid turbulence. In Helium, 
the formation of a vortex tangle is assumed to follow the 
onset of an instability in the vortex array 0]. It has re- 
cently been suggested that an analogue of this so-called 
Donnelly-Glaberson instability may be relevant for neu- 
tron stars [5|, i6j] . If this is the case, one would ex- 
pect it to have interesting repercussions for neutron star 
precession. In this Letter we conffi'm this expectation 
by demonstrating that short-wavelength instabilities are 
generic in precessing superfluid neutron stars. 

Plane wave analysis. — Our main objective is to in- 
vestigate whether analogues of the Donnelly-Glaberson 



instability are likely to occur in a neutron star interior. 
Our analysis is based on the standard two-fluid picture, 
where the superfluid neutrons are dynamically distin- 
guished from a conglomerate of comoving superconduct- 
ing protons and normal electrons. We will loosely refer 
to the latter as the "protons" in the following. Variables 
associated with each fluid will be labelled by x = {n, p}. 
The smooth-averaged hydrodynamics of the system is 
governed by two coupled Euler-type equations, see Q for 
more details. In a frame rotating with angular velocity 
ri* we have 
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Here the fluid velocities are denoted by v]^, we have in- 
troduced the convective derivatives — dt + v^Vj and 
Xp = Pp/Pn is the density fraction. The scalars tp^ are 
the sums of specific chemical potentials and the gravita- 
tional potential 0]. For simplicity, we assume that both 
fluids are incompressible, i.e. we have V^z;* = 0. In the 
interest of clarity, we also ignore the entrainment effect in 
this study. A key property of the system is that neutrons 
and protons are coupled via mutual friction, a force 
mediating the interaction between the quantized neutron 
vortices and the proton fluid/magnetic fluxtubes. The 
standard expression for this force is, see 0], 
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where = vl^ — and the neutron vorticity is given 
by LjJ, — e''^''WjV^. A "hat" denotes a unit vector. This 
form for the mutual friction force results from balancing 
the Magnus force that acts on the neutron vortices and 
a resistive "drag" force which represents the interaction 
between the vortices and the charged fluid Repre- 
senting the drag force by a dimensionless coefficient TZ, 
one finds that 
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In the most commonly considered case, the mutual fric- 
tion arises from scattering of electrons off the vortex's 
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intrinsic magnetic field. This leads to a relatively weak 
couphng, with 0, 0| , 
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It is, however, not established that it is this limit that 
applies. Hence, one must also consider the case of strong 
coupling which follows from taking TZ ^ 1. This trans- 
lates into 
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The strong coupling limit is relevant if the interaction be- 
tween neutron vortices and fluxtubes is efficient [sl. [ol. \ld^ . 
if there is a fluxtube cluster associated with each neutron 
vortex 11 1, or if there is significant vortex pinning [l^ 
(in the limit TZ —>■ oo the vortices can be considered as 
perfectly "pinned" ) . 

Returning to the Euler equations ((iJ-Q, only the pro- 
ton equation contains a shear viscosity term. This is 
because the dominant process is expected to be electron- 
electron scattering. The upshot of this is that the neutron 
fluid is not directly affected by shear viscosity. The rel- 
evant viscosity coefficient has been estimated in flst . 
For a uniform density star with M — IAMq, i? = 10 km 

0.1 (the canonical values we will use later) we 
i 10^(T/108K)^2 (,j^2/g_ 

We consider perturbations of cqns. H]) and ([2]) for a 
background configuration where both fluids rotate rigidly 
with w^o = e'-'''(ri^ — ilj)xk- By allowing for an arbitrary 
orientation of the angular velocity vectors, this configu- 
ration can represent the standard free precession modes 
of a two-fluid star IJ, |l5|. We then linearise the Eu- 
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ler equations, focussing on short-wavelength motion by 
making the standard plane-wave decomposition 
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and similarly for all other variables. Since we expect the 
flow along the background vortex array to play a central 
role [4, 5], we carry out the perturbation calculation in 
the neutron frame. That is, we take H.^ = fl^^ = rinfi^. 
In order to simplify the problem, without any real loss of 
generality Q, we only consider waves propagating along 
the vortices, i.e. fc* — Then the fact that we have 

assumed the fluids to be incompressible means that the 
waves are transverse, fiiAl- = 0. After some algebra, cf. 
for a similar analysis, the perturbed versions of ([T]) and 
([2]) lead to a 4 x 4 system, the determinant of which pro- 
vides the dispersion relation for short-wavelength waves. 
A detailed analysis of the problem will be provided else- 
where. Here we focus on the modes that may become 
unstable. 

Let us flrst consider the weak drag limit. Then we flnd 
a mode with frequency (with viscous corrections of order 

lAii) 
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Here represents the relative linear flow along the 
(background) neutron vortex array. In our case we have 
■w\\ — — nVyfefi^a;*^. In the local analysis we have taken 
w\\ to be constant. Hence it is clear that our analysis 
is only consistent for short wavelength motion. Anyway, 
from ^ we see that the system is unstable (Im ct < 0) if 
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As discussed in Q, the solution ([8]) represents inertial 
waves in the neutron fluid. This instability is the exact 
analogue of the Donnelly-Glaberson instability in Helium 
0, and hence its existence should come as no real sur- 
prise. As in Helium, one would expect the onset of the 
instability to lead to the formation of tangled vortices, re- 
connection and superfluid turbulence. Since turbulence 
alters the form of the macroscopic mutual friction force 
[3, 01 , it is not yet clear how the system will evolve once 
the unstable waves grow to large amplitude. 

As far as we are aware, the strong drag problem has 
not been considered previously. Interestingly, there are 
unstable modes also in this case. The nature of the in- 
stability is, however, more complex. In the strong drag 
limit, with B — and B' = I, we flnd a mode with fre- 
quency 
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This result clearly shows that there will be unstable 
waves (representing coupled inertial waves in the neu- 
tron/proton fluids). In the inviscid {vcc = 0) limit the 
instability is active provided that 
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As in the weak drag case, one would expect the onset of 
this instabifity to lead to tangled vortices. 

Implications for precessing neutron stars. — In order 
to discuss the implications of the above results we need 
to make contact between our background configuration 
and the global precession motion. Fortunately, this is 
straightforward. The precession of a two-component neu- 
tron star model, including mutual friction coupling, has 
already been discussed in [3]. The simplest model con- 
sists of two components that rotate rigidly. The neutron 
component is assumed spherical with moment of inertia 
In- At the same time, the protons (including the crust) 
are assumed to be slightly deformed in such a way that 
Ip — h — h = I3 / + e) with e <§; 1 (in a principal co- 
ordinate system where the deformation axis is along £3). 
When perturbed away from alignment of the two rota- 
tion axes, r2^, the crust precesses with a certain frequency 
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and observable wobble angle 9w (the angle between the 
deformation axis, is, and the total angular momentum 
axis) [l5| . 

The plane-wave analysis is consistent for the precess- 
ing system provided that the two rotation vectors fij^ can 
be considered fixed. This is true as long as the preces- 
sion period Ppr is significantly longer than the timescale 
associated with the local waves. As already mentioned, 
the wavelength of the waves we consider must also be 
short enough that w\\ can be treated as a constant. If 
these conditions hold then we are simply considering lo- 
cal perturbations of a given precession model. In order 
the check whether this system is locally stable we only 
need to work out w\\ from the precession solution. If an 
instability is present, then the precession solution must 
be considered questionable. It certainly cannot be the 
case that the two components rotate rigidly, a key as- 
sumption in the standard analysis [l^ . 

Weak drag slow precession. — In the weak drag limit, 
there exists a long period precession solution that is 
slowly damped by mutual friction. For this solution we 
have 



precessor PSR B1828-11 [H the spin period is 0.4 s and 
in order to have precession with the observed period one 
would need e « 10^^. It is then clear from that 
precession with a wobble angle of a few degrees is safely 
in the stable regime. Nevertheless, the weak drag result is 
not without interest. Consider for example a millisecond 
pulsar with a maximally strained crust. From (|16[) we 
see that if the spin period is 1 ms, then precession with 
larger than a degree would be unstable provided that 
e > 6 X 10"''. Since the theoreticall y p redicted range for 
crustal deformations has e < 10~^ [16|], we see that our 
result puts a constraint on slow precession in very fast 
spinning neutron stars. 

Strong drag fast precession. — In the stron g dr ag limit, 
the relevant precession solution is such that [14| 



and we find that 



14| 



2tt9iuX2 In 

p z 



(17) 



(18) 



Ppr « P/e , and td 



27ri3e 



(12) 



where Ppr is the precession period, id is the damping time 
are P is the rotation period of the star. We then find that 
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where X2{< R) is one of the coordinates associated with 
crust system. This estimate can be used in ^ to show 
that all waves with wavelength (A = 2n/k\\) shorter than 
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are unstable. However, there must be a short wavelength 
cut-off for the instability. To make progress it would 
seem natural to assume that our analysis becomes invalid 
once the wavelength is so short that the fluid description 
is no longer relevant. Then it seems reasonable to use 
something like 
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where d^ is the intervortex spacing. Since we need to 
have Aniax > Amin in order to argue that the instability 
is relevant we see that we must have 
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What does this result tell us? It suggests that, if the 
drag is weak, the superfluid instability is unlikely to play 
a role in slowly spinning systems. For the archetypal 



The (inviscid) instability criterion (jll[) then implies that 
waves with wavelength shorter than 
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(we have assumed /p//n ~ a;p) will be unstable. 

However, as is clear from (ITU)) , the unstable strong drag 
modes are affected by viscosity. To unveil the detailed 
behaviour we have solved the dispersion relation numeri- 
cally for a range of parameter values. Typical results are 
shown in Figure [1] This figure shows Tg^ow as a function 
of A and illustrates how the importance of shear viscosity 
varies with temperature. The results for core tempera- 
tures T — 10^ K and T — 10^ K show a clear transition 
from a regime where the inviscid approximation to (|10p 
is valid (above A « 20 cm and 10^ cm, respectively) to a 
short wavelength regime where viscosity alters the result. 
However, a surprising feature appears as one proceeds to- 
wards shorter wavelengths for a fixed temperature. As fcy 
becomes large, it turns out that there is a cancellation 
of the leading order viscosity terms, cf. PH)) . For short 
wavelengths, the mode frequencies are instead accurately 
(with errors of order \/k\\) described by ([5]) with B' = 1. 
Hence, for A Amax the short wavelength modes grow 
on a timescale given by 
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FIG. 1: Growth timescales rgrow(-^) for unstable superlluid 
waves in a precessing neutron star. This particular model has 
TZ = 10'', representing the strong drag regime, P = 1 s and 
9m = 1°. The dotted horizontal line shows the fast preces- 
sion period Ppr that follows if we take Xp = Ip/In = 0.1. We 
compare our numerical results for three different core temper- 
atures (solid lines) to two approximations. The short wave- 
length approximation and the 7?. 3> 1 approximation (|10p are 
both indicated by dashed lines. The region where a short- 
wavelength instability operates (I) for T = 10^ K is high- 
lighted. Finally, the A < Amin region where we assume that 
the hydrodynamical description fails is shaded. 



well within the stable regime. In contrast, our results 
have serious implications for systems in the strong drag 
regime. We predict that these systems will suffer local 
instabilities, possibly leading the formation of superfluid 
turbulence, for a wide range of the relevant parameter 
space. This calls into question the standard precession 
model, which is based on two co-existing fluids rotating 
as solid bodies [IJ], and any conclusions drawn from it. 



In particular, one would note Link's argument that 
the coupling between vortices and fluxtubes ought to lead 
to fast precession according to (fT7| . Since this is contra- 
dicted by the observed slow precession of PSR B1828-11, 
Link suggests that our understanding of the neutron star 
core physics is wrong and that the protons would actu- 
ally form a type I superconductor (without fluxtubes). 
Our results add an element of doubt. We have essen- 
tially shown that the strong drag fast precession solution 
may be inconsistent for a neutron star spinning at the 
rate of PSR B1828-11. If the precessing motion trig- 
gers a range of unstable short wavelength waves then the 
original solid-body assumption that led to (fT7|) . cf. [3], 
cannot hold. The precession problem may thus be more 
complex than usually assumed, and a consistent descrip- 
tion of fast precession must properly include superfluid 
wave dynamics and potential turbulence. 

This work was supported by PPARC/STFC via grant 
number PP/E001025/1. 



For consistency the unstable waves need to grow on a 
timescale that is short compared to the precession period. 
If we require (say) Tgrow < O.lPpr, then we have 

The corresponding instability region is indicated by an I 
in Fig. [T] Moreover, in order to have A > Amin (noting 
that the short wavelength cut-off remains as in the weak 
drag case) we must have 

This shows that the short-wavelength instability con- 
strains a wide range of fast precession models. From 
the results in Fig. [1] it is also clear that there may exist 
a medium wavelength instability regime (well approxi- 
mated by (fTOp . This is relevant for temperatures above 
10^ K, and could well lead to the fastest growing insta- 
bility in young neutron stars. 

Brief discussion. — In this Letter we have demon- 
strated that short wavelength superfluid instabilities may 
operate in freely precessing neutron stars. In the weak 
drag regime, the instability affects only rapidly spin- 
ning stars that have significantly deformed crusts. PSR 
B1828-11, the currently best candidate precessor, lies 
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